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Abstract. This paper addresses the Paris barrier options of [CJY] and 
their valuation using the Laplace transform approach. The notion of the 
Paris barrier options is extended such that their valuation is possible at any 
point during their lifespan, and the pertinent Laplace transforms of [CJY] 
are modified when necessary. 
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1. Introduction: An interesting application of the theory of Brownian excursion the- 
ory has been observed in [CJY] in connection with so-called barrier options. These are 
widely traded options which come as puts or calls that are activated or deactivated as 
soon as their underlying hits a prespecified barrier level. As explained in [CJY], this 
abruptness entails a number of undesirable consequences and problems both on a theo- 
retical and practical level. The idea of the Paris barrier options of [CJY] is to cushion 
this abruptness. This is by introducing a systematic delay for the consequences of hitting 
the barrier to become effective. And the insight of [CJY] is that well adapted means for 
this are afforded by excursion theory in general and the Brownian meander and Azema's 
martingale in particular. 

Still, Paris barrier options by now seem to be traded by the leading institutions only and 
not in large quantities. And the author is grateful to Pliska (UIC, Chicago) for pointing 
out to him a principal difficulty with the concepts of [CJY] . Trading Paris barrier options 
will typically lead to situations where the underlying has been staying above or below the 
respective barrier for a certain period of time but where this excursion has not yet lasted 
long enough for action to be taken. And there seem to be no results in [CJY] how to value 
Paris barrier options in such situations. In fact, such situations do not seem to be covered 
by the very constructions in [CJY]. One aim of this note so is to provide the necessary 
conceptual extensions to meet this criticism. 

For valuing Paris options, Laplace transforms of the Paris option densities have been com- 
puted in [CJY]. Pursuing Pliska's suggestions we moreover found it necessary to modify 
one of these two Laplace transforms. So I am very grateful to Pliska for his remark. And I 
also wish to thank Yor for the interest he showed in this work and for making available to 
me corrections to [CJY] collected by Jeanblanc-Picque which have also been incorporated 
in the exposition of this paper. 
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2. Black— Scholes framework: The analysis of this paper is in the Black-Scholes 
framework and uses the risk-neutral approach to valuating contingent claims. In this set- 
up there two securities. First there is a riskless security, a bond, that has the continuously 
compounding positive interest rate r. Then there is a risky security whose price process S 
is modelled as follows. Start with a complete probability space equipped with the standard 
filtration of a standard Brownian motion on it that has the time set [0, oo). On this filtered 
space one has the risk neutral measure Q, a probability measure equivalent to the given 
one, and a standard Q-Brownian motion B such that S is the strong solution of the 
following stochastic differential equation: 

dS t = (r-5)- Sfdt + a- Sf dB t , t e [0, oo). 

Herein the positive constant a is the volatility of S. The constant 5 depends on the security 
modelled. For instance, it is the dividend rate if S is a stock. 

3. Paris options: This section takes up the Paris barrier options of [CJY, §2]. Their 
new idea for cushioning the impact of the underlying hitting the barrier is as follows. They 
require the underlying S to spend a minimum time D > above or below their prespecified 
barrier L > before the option is knocked in or knocked out. 

The Paris down-and-in call to be considered in the sequel is the following European-style 
contingent claim on S written at time to and with time to maturity T. Its payoff at T is 
that of a call on 5 with exercise price K: 

{S T ~ K)+ = max{5 T - K, 0} , 

if S t is less than L during a connected subperiod of length at least D of the monitoring 
period [to, T]. Equivalently, there is a point in time a such that the interval I = (a,a + D) 
is contained in the lifespan [to, T] of the option and S t < L for any t in I. Otherwise the 
call expires worthless. 

Now let t be any time in [to, ^] at which the Paris down-and-in call has not yet been 
knocked in. The time-t value Cd,i of this call is then a function of K, r = T—t, So, r, 5, L 
and .D, and mostly any of these arguments is suppressed from the notation. The valuation 
problem depends on an excursion time H^ t defined in the sequel such that the Paris 
down-and-in call is knocked in before its maturity T if and only if H^ t < T. Granting 
this, using the arbitrage pricing principle, the time-t price Cd,i of the Paris down-and-in 
call is then given by the conditional (J-expectation: 



Cd,i = e rT E t 



<j}{S T ) ■ l{JT iit <T}] ■ 



with 4>(St) equal to the call payoff (St — K) + , and with 1{h l t <T} the characteristic 
function of the event H L t < T. 

In defining H Lyt at time t two constellations are to be distinguished. Namely, either today's 
security price St is equal to or above the barrier L, or it is below L. In the first case with 
St > L, the Paris down-and-in call is knocked in before its maturity T iff S remains 
smaller than L for all points in time of an interval of length at least D contained in [t, T]. 
Thus let H L t in this case denote the first time s greater than or equal to t + D such that 
S u < L for any time u in the interval (s — D, s). 
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In the second case where St < L, the security price S has already been staying below L 
for some connected period time during the lifetime of the option. The down-and-in call is 
then knocked in before its time to maturity T if the following happens. The security price 
S continues to stay below L also during the period of time from today until the point in 
time t + d with d smaller than D, i.e., one has S u < L for all u in (t, t+d). Thus define 
Hl i = t + d in this case. If the level L is hit by S before time t+d, however, the clock 
for the minimum lenght D is restarted. So define H^ t in this case to be the first time s 
greater than or equal to t+D such that S u < L for any u in the subinterval (s — D,s) of 
the positive real line. 

The discussion of Paris barrier options would not be complete without having referred to 
[CJY, §2] for details on the whole family of Paris barrier options. Calls and puts with 
either of the following barrier types: down-and-out, down-and-in, up-and-out, up-and- 
in. Valuing these eight types of options is reduced in a standard way to the case of the 
above Paris down-and-in call. Indeed, put-call parities are given in [CJY, §5], and [CJY, 
§4.2] reduces the valuation of the out-calls to that of the in-calls, mutatis mutandis. This 
paper thus concentrates on valuing the Paris down-and-in call, in the sequel also referred 
to as the Paris option for simplicity. 

4. The basic valuation identity: Valuing the Paris option in the Black-Scholes con- 
text is made precise by the following basic valuation identity for its time-t price Cd,i 



for any reals u > and y, where w = a -1 • (r— 5— a 2 /2), with (3(S t ) = a -1 log(K / S t ) , where 
— Hi^ — t, and with //* the measure for Brownian motion at time H^. The expectation 
is taken with respect to a Girsanov transformed risk neutral measure Q* constructed in §6 
when proving the basic valuation identity. Moreover notice that by construction, hb(u, ) 
is zero for u less than D, respectively d, depending on today's situation. 

5. Results on the Laplace transforms of the densities: This section discussed he 
Laplace transforms with respect to time of the Paris option density hb of §4. For any fixed 
real number y, recall they are defined by: 



for any complex number z with sufficiently big real part. Moreover define the function \1/ by: 




with the Paris option density hb given by: 






4 



Michael Schroder 



for any complex number z, and choose on the complex plane with the non-positive real 
axis deleted the square root defined using the principal branch of the logarithm. The 
computations of [CJY, §§5, 8] then essentially give the following two results. 

Proposition A: For any real number y, the Laplace transform of hb( , y) is a holomor- 
phic function on the right complex half-plane. 

Proposition B: Suppose b = a -1 log(L/S t ) is non-positive. For any fixed real number 
y, the Laplace transform of hb( , y) is given by: 



b 



V2D 



for any complex number z with positive real part 



^_ f°° . -£-\b-x- V \yfc 

D \/2Dz ty(V2Dz ) 



If b = a~ l \og{L/S t ) is positive, i.e., when today's security price S t is below the barrier L, 
we have modified the Paris option of [CJY]. Today's situation is then such that the price 
of the security has already been staying below the barrier L for a certain connected period 
of time during the lifetime of the Paris options. For the Paris option to be knocked in 
the price thus has to stay below the barrier L only for another connected period time of a 
length d < D from today, i.e., S u < L for all u in [t, t+d). In this way, Paris options now can 
be valued at any point of their monitoring periods. The desirability of this modification 
has been pointed out to me by Pliska. I regard the following modification of the valuation 
results of [CJY] as a direct result of his suggestions. 

Proposition C: Suppose b = a~ 1 log(L/St) is positive. For any real number y, the 
Laplace transform of hb( ,y) is given on the right half-plane as the following four-term 
sum of Laplace transforms: 

S*(M ,y)) =Erfc (-|=) • ^(fe M ( ,„)) + -J= . £'(h K2 ( ,»)) 

+ Erfc (db) Erf fob) ■ ^ { • v) ) + Erf (db) ' • ^ 



2d' \V2ds V / VV2d 

Herein h^k are the functions on the positive real line times the real line defined by: 

(y-b) 2 



u ( \ i ( \ 1 f Vu-d 2 (u-d) 

hb ^ y) = 1 ^ {u) -^dTD\-^- e 



_ 1 (y-b)VD c -2ife&) Erf / (y-b)VD \ 1 
2 ^/u-d+D \J2(u-d)(u-d+D) J J 



\^2(u-d)(u-d+D) 

jl_ (y-2bf 

h b ,±{u,y) = l {d)0o) {u)-^== ■ f b ,u(y) where fb,u{v) = e 



Tu, g 2u 
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for any positive real number u and any real number y, while the remaining two functions 
are defined using their Laplace transforms as follows: 

nd ~ zw / poo _ x 2 _ \ 

&{{hb,2{ ,y))(z) = / -=—==-( e x v 2z f bid (x)dx)fi b (dw), 
Jo \ / 2zW(v2Dz) \Jr / 

for any positive real number y and any complex number z with positive real part. Herein 
Hh denotes the law of the first passage time to the level 6, given on the positive real line by 



Hb(dw) = 




6. Proof of the basic valuation identity: For proving the basic valuation identity, 
consider the restarted at time t Brownian motion B*(u) = B(t-Hi)—B(t) and the normalized 
Brownian motion with drift W* given by 



w: = -io g ( 

a \ 



S-, 



t+u 



= -oju + B* , u e [t, oo) . 



By construction we then have S t+U < L iff W*(u) < b and H^ t < T iff H£ < r = T — t 
Using Strong Markov, thus transcribe the price Cd,i of §3 of the Paris option as follows: 



S t .e W ^-K 



) ■ 1 {h;<t}_ 



Now proceed in analogy to [CJY, §§4,5]. Using the particular case [KS, p.l96f] of the 
Cameron-Martin-Girsanov theorem, change measure from Q to Q* such that W* becomes 
a Q*-Brownian motion. Iterating the Q*-expectation thus obtained for Cd,i, we have 



Cd,i 



E* 



1 {H *< t} E* 



(St-e 



<jw: 



-K 



+ 



Apply the strong Markov property of Brownian motion using as ^"-stopping time. 
The resolvent so obtained further simplifies since the random variables and W*(H£) 
are independent. With /j,* the measure for W*(H£), the basic valuation identity follows. 



7. An intermediate identity for the proofs of the Laplace transforms: 

the concepts introduced in §4, we have the following intermediate identity: 



With 



Lemma: For any real number y, one has 

&(h b ( ,y))(z) 



[ E* 


r 


e 


/R 





\x-y\V2z 



fi*(dx) 



for any complex number z in the right complex half-plane. 
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For proving the Lemma interchange the Laplace transform with the two exponential inte- 
grals defining the Paris option density hb(u, y) in §4 to get: 



3(h b ( 7 y))(z)= / E* 



1 



HT y/2*{u-Ht) 



(x-y) 2 
■e 2 ( u ~ H t) du 



/j,*(dx) . 



In this last integral successively change variables w = u—H£ and separate the expectation 
from the Laplace transform to get: 



L 



E* 



-zHT 



V2 



e )(z)fi*(dx). 



nu 



Using e.g. [D, Beispiel 8, p.50f] for the Laplace transform then completes the proof. 

8. Computing the Laplace transform in the case b equal to zero: The compu- 
tation of the Laplace transforms of §5 for /i b is by reduction to the case 6 = where we 
have to show the key relation: 



E* 



^(zVD) = tf(0) = 1, 



for any complex number z with positive real part. This section reviews the key ideas of its 
proof in [CJY, §8]. The computations are based on the Azema martingale and properties 
of the Brownian meander. The Azema martingale is the martingale on [0, oo) for the slow 
Brownian filtration F + given for any t > by: 

IH = (sgnWt)Vt-gt- 

Herein g t is defined as the supremum over all real numbers s < t such that W*(s) = 0. 
The process given by: 



m t (u) = 



1 



y/t-gt 



\W*(g t + u(t-g t ))\, ue [0,1] 



is a Brownian meander and is independent of the a-subfield J r+ (gt) of the slow Brownian 
filtration. Its law is independent of t and given by xexp (— x 2 /2) ■ l(o,oo) dx- The tautology 
W*(t) = m t (l) ■ [it thus implies the following identity: 



E* 



zWT- 



for any positive real number z. Fixing any such z, the crucial point is to exhibit the product 
^{—z-h{Hq)) exp (—z 2 Hq/2) as a martingale. This is achieved by a boundedness argument 
showing the argument of the above conditional expectation to be uniformly integrable for 
any time t up to Hq . Using another form of the optional stopping theorem it follows that 
the Q*-expectation of the martingale stopped at Hq is equal to the Q*-expectation of it at 
time zero, whence to \I/ (0) which equals one. With the random variables Hq and W*(Hq) 
independent, the key relation results. This is an identity between functions holomorphic 
in particular on the right half-plane. Using the identity theorem, it thus remains valid for 
any complex number z with positive real part. 
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9. Computing the Laplace transform for b non— positive: Using §7 Lemma the 
proof of §5 Proposition B reduces to compute the expectation of exp (— zH b ), for any 
complex z with Re (z) > 0, and the density /i* . Following [CJY, §8.3.3], this is by reduction 
to the results of §8 using restarting arguments. Indeed, decompose H£ as follows: 

H£=T b + H**. 

Herein T b is the first passage time of W* to the level 6, and Hq* is defined as follows. It 
is the smallest non-negative point in time s at which the restarted-at-time-T^ Brownian 
motion W**(u) = W* (T b +u) — b is zero for the first time after having been less than zero 
for a connected period of time of length at least D. 

To compute the expectation required, notice that T b and Hq* are independent random 
variables. The conditional expectation at time T b of exp ( — (w 2 /2)H b *) thus is the prod- 
uct of exp (-{w 2 /2)T b ) and the expectation of exp (-{w 2 /2)Hq*). With the law of Hq* 
equivalent to that of Hq, this last expectation is given by the key relation of §8. So it is 
deterministic in particular. Take expectations to time zero of this product. The expecta- 
tion of exp (— (w 2 /2)T b ) that remains to be computed is the Laplace transform of the law 
of T b at w 2 /2. Using [D, Beispiel 8, p.50f] for instance, it is standardly seen to be equal 
to exp(-|6|(2 • (w 2 /2)) 1 / 2 ) and thus to exp (+&«;). With w = (2z) 1/2 one gets: 



E* 



l 



V(V2Dz ) 



E* 



-y/2z T b 



1 



by/2z 



ty(V2Dz ) 



To determine \i* notice the tautology W*{H*) = W**(Hq**) +b. Thus W*(H*) < x iff 
W**(Hq*) < x-b. From §8 the law of W**(Hq*) is minus that of mi, whence 



Q* {W* {Ht) e dx) = 1 ( _ 00)6] (x) -{b-x)-e 



2D 



dx 



Substitute into §7 Lemma and change variables w = b—x to obtain §5 Proposition B for any 
complex number z with sufficiently big positive real part. Extend to the right half-plane 
using the consequence of the key relation of §8 that \E , ((2/J^) 1 / 2 ) has no zeroes there. This 
completes the proof of §5 Proposition B and the 6-non-positive part of §5 Proposition A. 

10. Computing the Laplace transform for b positive: The proof of §5 Propo- 
sition C modifies the argument of [CJY, §8.3.3]. Using §7 Lemma, the first task is to 
compute the density n* and the Q*-expectation of exp(— zH£) for any complex number 
z with positive real part. The two cases where the level b is hit before or after time d are 
to be distinguished. Thus decompose the underlying probablity space Q into the set A on 
which the first passage time T b of W* to the level b is less than or equal to d 

A = {uj e O | T b (oj) < d} 

and its complement O \ A on which T b is bigger than d. This induces the decomposition: 



E* 



e~ zH » 



E* 



I A ■ e 



+ E* 



of the Q*-expectation of the random variable exp(— zH^). By construction H£ = d on 
the complement of A. So the second above summand is exp (—zd)Q*(T b > d). On the set 
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A the level b is reached before the critical time d and the clock for the excursion is reset 
to zero at the first passage time T b to the level b. Accordingly one has the decomposition: 

IT* rp | TT** 

H b = l b + li Q , 

with Hq* defined as follows. It is the smallest s > at which the restarted-at-time-Tb 
Brownian motion W**(u) = W* (T b -Hi)— b is zero for the first time after having been less than 
zero for a connected period of time of lenght at least D. With this random variable Hq* 
independent from T b , the conditional expectation of the random variable 1a- exp (—zH£) at 
time T b thus equals exp (— zT b ) times the expectation of exp (—zHq*). The key relation of 
§8 now applies to Hq* and identifies this last expectation as the reciprocal of ^(^Z-D) 1 / 2 ). 
The expectation of 1a • exp (— zT b ) on the other hand, is given by integrating exp (— zw) 
from zero to d against the density fi b of T b . Summarizing, it so follows for any complex 
number z with positive real part: 

E* [ e -*^] = Q*(T b > d)e~* d + ^^=y • 

For determining ji* decompose the random variable W*(H b ) with respect to A: 

W*(H* b ) = (W**(H**)+b)-l A + W*(d)-l Q \ A - 

Using the independence of T b and H$*, the law /j* A (dx) = Q* ((W** (Ho**)+b) G dx; T b < d) 
of the first summand is obtained as the convolution of the laws of T b and W**(Hq*): 

£ dx 



Mdx) = Q*(T b < d) ■ l ( _ 00)&] (x) • (b-x) ■ e ™ [} 

For the law fi^ A (dx) = Q* (W*(d) G dx; T b > d) of the second summand notice T b > d iff 
W* (t) < b for all t < d. Thus the first passage time condition T b > d transscribes into the 
condition that the running maximum max{VF*(t)|t < d} is smaller than b. Using [H, p. 9], 

Applying §5 Lemma, the Laplace transform i? > (/ife( , y))(z) of h b ( ,y) at z is equal to: 

<?(h b ( , y)){z) = E* [e-^] 6 lX Jf^ (»* A + V&\ A )(dx) 

using that the expectation factor is deterministic and independent of the variable x. So it 
suffices to compute the two integrals that make up the second factor of this product. For 
the first of these, on substituting for /i A and changing variables w = b—x, one gets: 

f e^-y^ Q*(T b <d) f°° e-^-^g 
/ == — u A (dx) = — / x ^= e ZL> dx . 
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For the second of these one analogously obtains 

\x-y\V2l , , -\ X -y\y/2Z x 2 (x-2b) 



2 

dx . 



C P ~\x-y\v ^ 1 r p -\x-y\v±z yx--iby 

Multiplying out, J/ ? (h b ( ,y))(z) is the sum of four terms of which each is seen to be a 
Laplace transform itself. Laplace inversion of this identity so is term by term. To identify 
the remaining functions h b ,3 and h b ,4, recall for instance using [D, Beispiel 8, p.50f]: 

/ 1 / — \ 1 -a 2 

JS*- 1 [ ^ • e~ av ^ ) (u) = -7=-e ^ 



valid for any positive real numbers a and u. The Laplace inverse of the improper integral 
factor of the above \i* A integral is, on applying Fubini's theorem, equal to: 



xe 2D ■ jS"" 1 (e- zd = \u)dx 



2z 

Successively apply the shifting theorem for the Laplace transform to take care of the factor 
exp (—zd) and the above Laplace inversion formula with a equal to \b— x — y\ to arrive at: 

1 1 r°° ( x 2 1 ( b - x -y) 2 \ 

.fn.\ / ~. e y2T) + 2(u-d) J d x 



1 If 



X 



This integral is seen to be the value of the function ^5 3 at u and y. Analogously, the 
Laplace inverse at any u > of the integral factor of the above f^^ A integral is equal to: 



■ p (x-y) 2 x 2 (x-2bY 

l( doo) (u)—= / e 2 M (e 2d -e 2d Wr . 



V2ivd ' ^J2ix{u—d) Jn 

Again suppressing the details of the computation, it is seen to be equal to the value of the 
function hb,4 at u and y. This completes the proof of §5 Proposition C. 
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